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It is shown that the infiiiite dimensional critical surface of general eucUdean lattice actions in a generic four- 
dimensional scalar field theory with interactions has a domain of special multicritical points where higher 
dimensional operators play a special role. RenormaUzed trajectories of higher derivative continuimi held the- 
ories with nontrivial interactions are traced back to special ultraviolet stable fixed points on the manifold of 
multicritical points. These fixed points have an infinite number of relevant directions which identify the uni- 
versaUty classes of critical higher derivative field theories. The relevance of the new fixed point structure is 
discussed within the context of the triviaUty Higgs mass bound. 



1. Introduction 

It is generally accepted that the only fixed 
point on the infinite dimensional critical sur- 
face of a scalar lattice field theory with in- 
teraction in four dimensions is an infrared sta- 
ble Gaussian one. In this scenario, the only 
renormalized trajectory in the codimension of 
the critical surface, which in Whson's termi- 
nology [1] defines a massive continuum scalar 
field theory, originates from the Gaussian fixed 
point when traced back to the critical sur- 
face, thereby describing a noninteracting sys- 
tem. Since the Gaussian fixed point is infrared 
stable, the massless scalar theory on the crit- 
ical surface is also trivial. This picture would 
suggest that the critical surface of the the- 
ory with internal 0(N) symmetry is a man- 
ifold of simple critical points separating an 
ordered (ferromagnetic, or broken symmetry) 
phase from a disordered (symmetric) phase. A 
trivial renormalized trajectory runs into both 
phases from the Gaussian fixed point. 

However, a scalar field theory with inter- 
actions and higher derivatives in the quadratic 
part of the Lagrangian creates a seemingly 
paradoxical situation. It is a continuum field 
theory in four dimensions with nontrivial in- 
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teractions [2]. In fact, a higher derivative term 
jij$n3<|)^ with a new mass scale M in the 
Lagrangian, renders the theory free of diver- 
gences. Since finite fluctuations cannot com- 
pletely screen the <!>'* coupling, renormalized 
interactions will not vanish in the continuum. 
These theories cannot be described by the triv- 
ial renormalized trajectories of the infrared 
stable Gaussian fixed point. 

One might suggest that M now plays the role 
of the cut-off, therefore nothing has changed 
in the triviality scenario. This, however, cannot 
be a satisfactory resolution of the puzzle. It 
is legitimate to investigate the critical proper- 
ties of higher derivative field theories when M 
is kept on the scale of the low energy spectrum 
and carmot directly play the role of the cutoff. 
This happens when M is kept fixed in inverse 
correlation length units of ordinary massive 
particles as the critical surface is approached 
(limg^cx; = const) and M becomes part 

of the low energy spectrum. By pushing the 
ratio higher and higher, M could ulti- 

mately take over the role of a regulator under 
some special circumstances. However, one can 
always make the fluctuations divergent again 
by adding new interaction terms, and the puz- 
zle will reappear. 

An explanation of the paradox will be of- 
fered by showing that the critical surface of 
a generic scalar lattice theory is more com- 
plex than suggested by the triviality scenario. 
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There exists a domain of special multicriti- 
cal points where higher dimensional operators 
play a special role. The connection with the 
triviality picture which dominates the largest 
part of the critical surface can be explained by 
a complex crossover mechanism. The concept 
of renormalization group flows will guide our 
discussions. 

2. Renormalization Group and Triviality 

We will apply Wilson's renormalization 
group analysis [1] to general effective Hanul- 
tonians of scalar field theories with a physical 
cutoff at a very large energy scale A. The the- 
ory at energies above A could be a different 
field theory, or something completely new. 
The focus of owr interest is the low energy de- 
scription, far below the cutoff scale. Below the 
cutoff A, the most general effective Hamilto- 
nian (euclidean action) in d dimensions is de- 
fined in terms of an infinite number of momen- 
tum dependent couplings w**^' {qi,..., qn), 

n=l J 

X<5<'*'(Si4i)*(4i)---*(4n). (1) 

The only constraints on the Harmltonian come 
from the symmetries of the system. For ex- 
ample, only even terms are allowed in Eq. (1) 
under the $ ^ Ising like symmetry. Al- 
though the implementation of the physical cut- 
off is quite arbitrary, we will assume a smooth 
momentum cutoff A in continuous euclidean 
space, or a lattice discretization at small dis- 
tances where the lattice spacing a is related to 
the lattice momentum cutoff by A = n/a. 

A renormalization group transformation in- 
volves the lowering of the cutoff from A to 
Ap and a rescaling of the field $, as re- 
quired by the critical properties of the sys- 
tem. The transformed effective Hamiltonian 
J-Cp{^) is described by a set of new couplings 

Wp"' . . . , ifn) which satisfy a coupled sys- 
tem of renormalization group equations. At a 
fixed point 5f * the Hanultonian is not chang- 
ing, and the low energy theory becomes inde- 
pendent of the cutoff. The renormalization 
group flow of the effective Hamiltonian 3-Cp 
around the fixed point is determined from a 



linear stability analysis. With 3-[p = J-f*+A3-[p 
we write 



P^AMp =L*{AMp) , 
dp 



(2) 



where I* is a linear operator. The solution of 
the linear renormaUzation group equations is 
given by 



yip 
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The eigenoperators O* of L*, with eigenvalues 
Aj = e^s are classified according to their scal- 
ing exponents yi, 

yi > relevant eigenoperator O* , 
yi = marginal eigenoperator Of , 
yi < irrelevant eigenoperator O* . 

Far away from the critical point, the expan- 
sion A3^p = I.iiJiip)0* is still applicable in 
terms of the complete set of eigenoperators 
Of. However, the p-dependence of the scal- 
ing fields pdp) will be governed by nonlin- 
ear renormalization group equations, replac- 
ing the simple exponential evolution. 

A Hamiltonian J-[ at p = is critical, if it 
converges to a fixed point with infinite correla- 
tion length, limp^oo.;^p = The couplings 
of critical Hamiltonians define the critical sur- 
face. The type of the critical behavior depends 
on the number of symmetry conserving rele- 
vant operators. For a normal critical point one 
has one relevant symmetry conserving opera- 
tor which determines the critical temperature 
in statistical physics, or the critical mass pa- 
rameter in euclidean field theory. At a mul- 
ticritical point one has several relevant sym- 
metry conserving operators, and consequently, 
several conditions for criticality. 

In order to explore the critical domain of the 
Gaussian fixed point, we start from a simple 
critical Hamiltonian, 



u^^Hq) = q^ + 0{q^) , 



(4) 



which flows to the Gaussian fixed point 

u^*\q) = q^in the limit of a sharp momentum 
cutoff implementation of the renormalization 
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group (other RG transfromations imply more 
complicated fixed point functions.) With scal- 
ing dimension d$ = ^ (rf - 2) for the scalar field 

i>, the eigenfunctions u% (qi) are homoge- 
neous polynomials of degree 2k in the momen- 
tum variables qi. The scaling exponents are 
given by yn,k = d - j(d - 2) - 2k. At d=4, the 
only relevant eigenoperator j d'^q^(q)^{-q) 
has the scaling exponent yzfl = 2 and it de- 
scribes the mass term. The only nomedundant 
marginal eigenoperator is given by J <I>^ {x)d'^x 
with scaling exponent 3/4,0 = 0. From the 

quadratic equation p^jL/4,o(p) ~ ^\fl(p) of 
the marginal direction, the solution ^4,0 (p) ~ 
1 /Inp flows with slow logarithmic rate to the 
Gaussian fixed point. All other eigenoperators 
are irrelevant. 

The triviality scenario assumes that no other 
fixed points exist on the critical surface sug- 
gesting a manifold of simple critical points 
separating the ordered phase from the disor- 
dered phase. Earlier analytic and numerical 
studies had been consistent with this picture 
[3-7]. We will show, however, that there ex- 
ists a region of special multicritical points on 
the critical surface which describes a new class 
of nontrivial continuum euclidean field theo- 
ries. This will be demonstrated by a block spin 
renormalization group analysis on a simple hy- 
percubic lattice in d dimensions. 

3. Fixed Points in the Multicritical Region 

Consider first the quadratic lattice Hamilto- 
nian 



(5) 



where a scalar field ^{ii) is associated with 
each lattice point n. The general hopping term 
K(f) is not restricted to nearest neighbor inter- 
actions. A large finite volume, with V = I^, and 
periodic boundary conditions <l>(n -f Lcyn) = 
<^(n), K(f + Lem) = K{f) are assumed in the 
calculation. The block spin transformation, as 
schematically depicted in Fig. 1, defines new 
fields <i>'(n') on blocked sites n' where each 
block has 2^* sites of the original lattice. The 
transformation function of the block transfor- 



Figure 1. The blocked fields <!>' are defined on 
blocked sites n' which are designated by open 
circles. 

mation is defined by [8] 

Ta,i;[*',*] = -|aX(*'(n')-^ S <^>(/))' -(6) 

n' Jen' 

The parameter a in Eq. (6) should not be con- 
fused with the lattice spacing which is set to 
one, so that the dimensional quantities are 
all measured in lattice spacing units. After 
k blocking steps the partition fvmction Z = 
remains unchanged, 



const 



J 



g'rj,^j,'[4.',4.]-5f(*) _ 



(7) 



where 'J'j^^'[<J>',<J>] is obtained from Eq. (6) by 
the replacements a ^ Uk = and 

b bk = h^. The transformed Hamiltonian 
5f(fe)(<j)') is calculated by integrating out the 
variables <l>(n) of the original lattice. 



const 



(8) 



where a general quadratic form 

= I^IJ^'*^' (a ')<!>' (4')<J>'(-4')(9) 

is preserved. In Eq. (9) I' = 1/2*^, and the 
Fourier transform!^ (4) = Y.f^~^^^'<(f) is used 
in each blocking step. 
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With the choice b = 2 '2 + ^), the initial criti- 
cal Hamiltonian with if (4) = z{q^)^ + ■ ■ ■ will 
converge to a fixed point, limfc^<x, 5f — 
linife^o, if(fc) iq) - K* (q), where 

In Eq. (10) a* = a{l - and K*{q) rep- 

resents a line of fixed points parametrized by 
a. Higher powers {q'^)^,... (irrelevant opera- 
tors) in the initial couplings K{q) oi the critical 
Hamiltonian have no effects on the fixed point 
function K* (q). The normalization constant z 
will have an important interpretation. 

When K*(q) is expanded in powers of q^, 
the fixed point Hamiltonian in four dimesions 
has the simple form 

5f* ~ j(z(^')3 + ■ ■ ■)^{qm-q)d^q (11) 

in the infinite volume limit, with eigenvalues 
A2,fe = 2^"^*^ for the quadratic eigenoperators 



) {q'r''^^^{qm-q) . (12) 



A very similar picture is reproduced by using 
a diff erential renormalization group method in 
the continuum, if by analogy with the scale f ac- 
tor f7=2"'2+^), a new scaling dimension = 
j{d - 6) is chosen for the field $. The scaling 
exponents, given by yn,k = d - j{d - 6) - 2k, 
identify three relevant quadratic operators at 
n = 2, 



^2,0 



^^{qm-q)d'^q, 
\^q'^^{q)<^(-q)d'^q 



y2,o = 6 , 

3^2,1 = 4 , 



iq'y^{q)^{-q)d''q, 3^2,2 = 2 



One of the relevant quartic eigenoperators at 
n = 4 describes the A j ^'^{xjd'^x interaction 
term in a general Hamiltonian, 



<^'^(x)d'^x , 3/4,0 = 8 



(13) 



In fact, y4,k is positive for k = 0, 1, 2, 3. For 
general n, there exist an infinite number of rel- 
evant eigenoperators with positive scaling ex- 
ponents. Their significance can be identified as 



new interactions which remain renormalizable 
by power counting in the presence of higher 
derivative propagators. 

The extension of the above analysis to 0(N) 
internal symmetry is straightforward. It is also 
important to note that an infinite class of new 
fixed points K*(q) = z(q^)'^ + ■ ■ ■ can be gen- 
erated with positive integers a > 2, if b = 

2^(2+°') and new scaling dimension d$ = ■|(d- 
2a) are chosen. Our most studied model has 
a = 3 which serves as an example in the dis- 
cussion. 

The newly foimd fixed points, as described 
by K*iq), are generalized ultraviolet stable 
Gaussian fixed points in field theory terminol- 
ogy, with an infinite number of relevant direc- 
tions. We will show now that the renormal- 
ized trajectories of higher derivative contin- 
uum field theories with nontrivial interactions 
are traced back to these special ultraviolet sta- 
ble fixed points on the manifold of multicritical 
points. 

4. Continuum Higher Derivative Theory 

Consider the higher derivative Lagrangian 

£ = _i<j>^(x)(m2 + □ +M-4n3)$„(x) 

-Ao(*«(x)*„(x))' , (14) 

with internal 0(N) symmetry, where a summa- 
tion is understood over a = l,2,...,N. From 
a conventional viewpoint, the higher derivative 
term M"^ acts in the inverse euclidean prop- 
agator as a Pauli-Villars regulator with mass 
parameter M. It represents a minimal modi- 
fication of the continuum model when its eu- 
clidean path integral is rendered finite in four 
dimensions [2]. The euclidean propagator is 
equivalent to the sum of three simple poles 
where a complex conjugate pair of ghost poles 
is added to original massive particle. 

Since this model was discussed before [2], 
only the connection with the new fixed point 
structure will be presented. First, the scalar 
field will be given a new scaling dimension in 
the euclidean functional integral by the trans- 
formation qpa = m'^^a where m, identified as 
the inverse correlation length after renormal- 
ization, is measured in units of a large momen- 
tum cutoff A. With the notation z = m'*/M^, 
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the couplings of the eucHdean Hamiltonian, 
when expressed in terms of rescaled fields g?a, 

can be written as u^'^Hq) = + m'^q^ + z ■ 
(q^)^ and m***' = Aq ■ m^. 

In the limit m ^ 0, when z is kept fixed, 
the critical surface is approached in the univer- 
sality domain of the fixed point (q) which 
has three relevant quadratic operators in qpa- 
Therefore, a renormalized trajectory can be 
chosen which is traced back to 



z(4^)Va(4)<Pc 



(-q)d'^q 



(15) 



in the ultraviolet limit. It describes the con- 
tinuum field theory of Eq. (14) with nontriv- 
ial interaction in the infrared region. The 
interaction term of Eq. (13) scales now as 
Aom^(cpo(CPo,)^ along the renormalized trajec- 
tory. 

The nontrivial running coupling constant 
can be demonstrated by conventional renor- 
malization methods providing a more com- 
plete understanding of the critical behavior [9]. 

5. Nontrivial Running Coupling Constant 

The calculation of the scale dependent one- 
loop jS-function in the broken phase will illus- 
trate the nontrivial continuum limit of the the- 
ory. In addition to N-1 Goldstone excitations 
*o( = TTa, a = 1, 2, ... ,N - 1, there is also a 
massive Higgs excitation cr, *jv = v -i- cr, where 
V designates the renormalized vacuum expec- 
tation value of $Ar. A massive complex con- 
jugate ghost pair appears in all channels, as a 
consequence of the higher derivative term in 
Eq. (14). The following scale dependent pre- 
scriptions specify the renormalization condi- 
tions. 

1. The renormalization prescrip- 
tion on the one-particle irre- 
ducible Goldstone 2-point function, 
^r"^(p2)|p2=o = 1, sets the cor- 
rect normalization for the low en- 
ergy Goldstone modes. 

2- (47)^r^"(p')lp^=M2 =M-4 defines 
the renormalized ghost pole param- 
eter M. 

3. A scale dependent renormalized 
Higgs mass miij) is defined by 



raa(^2) = Zl^^ + Z3|J^/M'^ + 'm^i^) 

in the one-particle irreducible Higgs 
2-point function. Zi and Z3 desig- 
nate renormalization constants in 
the counterterms of the derivative 
parts of the euclidean Hamiltonian. 

4. The renormalization condition 5v = 
keeps the tree level relation v = 

2 

^ exact in every order. This de- 
fines a scale dependent renormal- 
ized coupling constant A(/i). 

The scale dependent one-loop jS-function is de- 
picted in Fig. 2. At low energies the jS-function 
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Figure 2. The one-loop jS-function is plotted 
against the logarithmic scale t = Ln^L/^/v^ from 
the numerical integration of the renormaliza- 
tion group equations. The point t = corre- 
sponds to iJ = V and t = on the logarith- 
mic scale. The initial condition A (0) = 10/24 is 
chosen with M/v = 14.3 which puts the ghost 
location into the multi-TeV range. 

is dominated by the Goldstone modes whose 
one-loop contribution is j^\^{t). Above the 
Higgs mass threshold the Higgs loop kicks in 
and the jS-fvmction becomes ^rA^ (t) which is 
the familiar one-loop form in the minimal mass 
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independent subtraction scheme. The onset 
of ghost effects becomes dramatic at the scale 
^ = where the ^-function turns around 
and begins to drop rapidly towards zero as the 
Goldstone and Higgs loop contributions are be- 
ing cancelled by ghost loops. The freeze of 
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Figure 4. The phase diagram of the lattice 
model at infinite bare coupling. The dotted 
line is calculated in the large-N expansion. The 
solid line displays the fixed z = m'^IM'^ ra- 
tio towards the continumn limit of the higher 
derivative theory. 



Figure 3. The rurming coupling constant 
is plotted from the numerical integration of 
the one-loop renormalization group equations 
with input parameters as defined earlier. 

the running coupling constant above the ghost 
scale is shown in Fig. 3 indicating that the sys- 
tem remains weakly coupled in the entire en- 
ergy range and the one-loop calculation is reli- 
able. 

6. Fate of the Lattice Higgs Mass Bound 

A hypercubic lattice structure was intro- 
duced [2,10,11] in non-perturbative com- 
puter simulations of the higher derivative La- 
grangian of Eq. (14). The phase diagram is 
shown in Fig. 4 in the limit of infinite bare 
coupling. Tuning the hopping parameter k 
to the critical line for fixed M corresponds 
to the triviality limit of the scalar field the- 
ory. In this limit, the dimension eight operator 
M'^^aU^^a becomes irrelevant and the dot- 
ted critical line with the exception of the end 



point at the origin represents a single univer- 
sality class. The continuum limit of the higher 
derivative theory without the underlying lat- 
tice structure is equivalent to the tuning of k 
towards the origin along a line of fixed z. A 
multicritical point is approached in this limit, 
as indicated by the sohd line of Fig. 4. The 
higher derivative term M"^<I>o,n^<I>o( becomes a 
relevant operator in this limit, and the theory 
is governed by the nontrivial fixed point. 

Higgs mass values from lattice simulations 
are shown in Fig. 5. In the 600 GeV to 700 
GeV range, the solid line of the simple hyper- 
cubic lattice action does not include any higher 
dimensional operators [3-5]. The dashed line 
corresponds to an F4 lattice action with dimen- 
sion six irrelevant operators [7]. A Symanzik 
improved simple cubic lattice action is also 
shown with dimension six irrelevant opera- 
tors towards the restoration of euchdean in- 
variance at finite correlation length [6]. In the 
much higher Higgs mass range, between 1.6 
TeV and 1.8 TeV, a dimension eight term is 
added as a relevant operator to the euclidean 
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Figvire 5. The Higgs mass bound of several lat- 
tice simulations is depicted. The arrows indi- 
cate one percent deviation from euclidean in- 
variance in a scattering amplitude as defined in 
[4]. 



action on a simple cubic lattice, to describe a 
strongly interacting Higgs sector with higher 
derivative Lagrangian [2,10,11]. The ghost lo- 
cation in this case is placed in the multi TeV 
mass range, and the lattice spacing could be 
completely eliminated. 

Large N calculations also indicate that at mn 
= 700 GeV the ratio M/niH is of the order of 
30 at infinite bare coupling. With the complex 
ghost location in the 20 TeV mass range, the M 
dependence in scattering amplitudes becomes 
practically invisible, well within the intrinsic 
ambiguity of the perturbative expansion. This 
is in sharp contrast with earlier lattice models 
[3-7] where the violation of euclidean invari- 
ance was larger than the ambiguity of the per- 



turbative expansion in the mn = 700 GeV mass 
range. We have to suggest, therefore, that the 
Higgs mass bound in the 700 GeV mass range 
in the earlier calculations was imposed by the 
artifacts of the underlying lattice structure. 

Since the complex conjugate ghost pair of 
the higher derivative Lagrangian theory evades 
easy experimental detection without violating 
unitarity [9,12], or Lorentz invariance, it serves 
as a model of the strongly interacting Higgs 
sector without technicolor, a scenario which 
was excluded in previous lattice studies. 
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